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Abstract —In this work we study the problem of rope 
sway dynamics control for elevator systems. We choose to 
actuate the system with a semi-active damper mounted 
on the top of the elevator car. We propose nonlinear 
controllers based on Lyapunov theory, to actuate the 
semi-active damper and stabilize the rope sway dynamics. 
We study the stability of the proposed controllers, and 
test their performances on a numerical example. 

I. Introduction 

Modern elevators installed in high-rise buildings are 
required to travel fast and ensure comfort and safety 
for the passengers. Unfortunately, the dimensions of 
such high-rise buildings make them more susceptible 
to the impact of bad weather conditions. Indeed, 
when an external disturbances, like wind gust or 
earthquake, hits a building it can lead to large rope 
sway amplitude within the elevator shaft. Large 
amplitudes of rope sway might lead to important 
damages to the equipments that arc installed in the 
elevator shaft and to the elevator shaft structure itself, 
without mentioning the potential danger caused for 
the elevator passengers. It is important then to be able 
to control the elevator system and damp out these 
ropes oscillations. However, due to cost constraints, 
it is preferable to be able to do so, with minimum 
actuation capabilities. Many investigations have been 
dedicated to the problem of modelling and control 
of elevator ropes [1], [2], [3], [4], [5], [6], [7], In 
[4], a scaled model for high-rise, high-speed elevators 
was developed. The model was used to analyze the 
influence of the car motion profiles on the lateral 
vibrations of the elevator cables. The author in [1] 
proposed a nonlinear modal feedback to drive an 
actuator pulling on one end of the rope. The control 
performance was investigated by numerical tests. In 
[6], a simple model of a cable attached to an actuator at 
its free end was used to investigate the stiffening effect 
of the control force on the cable. An off-line energy 
analysis was used to tune an open-loop sinusoidal 
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force applied to the cable. In [5], the authors proposed 
a novel idea to dissipate the transversal energy of 
an elevator rope. The authors used a passive damper 
attached between the car and the rope. Numerical 
analysis of the transverse motion average energy was 
conducted to find the optimal value of the damper 
coefficient. In [7], [8], [9], the present author studied 
the problem of an elevator system using a force 
actuator pulling on the ropes to add control tension 
to the ropes. The force actuator was controlled based 
on Lyapunov theory. Although, the sway damping 
performances obtained in [7], [8], [9] were satisfactory, 
we decided to investigate other forms of actuation and 
control. Indeed, in [7], [8], [9] we proposed a force 
control algorithm to modulate the ropes tensions, using 
an external force actuator, which was introduced at 
the bottom of the elevator shaft to pull on the ropes. 
However, retrofitting existing elevators with such a 
cumbersome actuator, can prove to be challenging 
and expensive. Instead, we propose to investigate here 
another actuation method, namely, using semi-active 
dampers mounted between the elevator car and the 
ropes. These actuators are less cumbersome than the 
typical force actuators and thus can be easier/cheaper 
to install. This type of actuation has already been 
proposed in [5]. The difference between the present 
work and [5], is that instead of using a static damper 
with constant damping coefficient tuned off-line, we 
use a semi-active damper, and use nonlinear control 
theory to design a feedback controller to compute 
online the desired time-varying damping coefficient 
that reduces the rope sway. We study the stability of 
the closed-loop dynamics, and show the performances 
of these controllers on a numerical example. One 
more noticeable difference with the work in [7], is 
that due to the semi-active actuation, the model of 
the actuated system is quite different from the one 
presented in [7]. Indeed, in [7] the model of the 
actuated system exhibited terms where the control 
variable was multiplied by the sway position variable 
only. In this work, the actuated model, exhibits terms 
where the control variable, i.e. the semi-active damper 
coefficient, is multiplied both by the sway position and 



the sway velocity variables, furthermore, the control 
variable term appears in the right hand side of the 
dynamical equations of the system, i.e. as part of an 
external disturbance on the system (refer to Section p|). 
These differences in the model, make the controller 
design and analysis more challenging than in [7]. 


The paper is organized as follows: In Section |If| 
we recall the model of the system when actuated 
with a semi-active damper mounted between the ropes 
and the elevator car. Next, in Section [TIIJ we present 
the main results of this work, namely, the nonlinear 
Lyapunov-based semi-active damper controllers, 
together with their stability analysis. Section IW] is 


dedicated to some numerical results. Finally, we 
conclude the paper with a brief summary of the results 
in Section lYl 

Throughout the paper, R, R + denotes the set of real, 
and the set of nonnegative real numbers, respectively. 
For x G M ;V we define x = y/x T x, and we denote 


by Aij, i = 1,.... n, j = 1, ...,m the elements of the 
matrix A. 


II. Elevator Rope Modelling 

In this section we first recall the infinite dimension 
model, i.e partial differential equation (PDE), of a 
moving hoist cable, with non-homogenous boundary 
conditions. Secondly, to be able to reduce the PDE 
model to an ODE model using a Galerkin reduction 
method, we introduce a change of variables and re¬ 
write the first PDE model in a new coordinates, where 
the new PDE model has zero boundary conditions. Let 
us first enumerate the assumptions under which our 
model is valid: 1) The elevator ropes are modelled 
within the framework of string theory, 2) The elevator 
cai' is modelled as a point mass, 3) The vibration in 
the second lateral direction is not included, 4) The 
suspension of the car against its guide rails is assumed 
to be rigid, 5) The mass of the semi-active damper is 
considered to be negligible compared to the elevator 
cai' mass. 

Under the previous assumption, following [3], [1] 
and [5], the general PDE model of an elevator rope, 
depicted on Figure [T] is given by 

P(w + + 2v ^WcTt + a lk) u hh t ) 

-^T(y,f)^^ + c p (| + n(f)^)n(i/,f) ( 1 ) 

+kdp{m + v( S)-§j J )u(y,t)8{y - l + l dp ) = 0 

where u(y, t ) is the lateral displacement of the rope, p 
is the mass of the rope per unit length. T is the tension 



Fig. 1. Schematic representation of an elevator shaft showing the 
different variables used in the model 


in the rope, which varies depending on which rope in 
the elevator system we are modelling, i.e. main rope, 
compensation rope, etc. c p is the damping coefficient 
of the rope per unit length, v = c -^~ is the elevator 
rope velocity, where l : R — » R is a function (at least 
C 2 ) modelling the time-varying rope length, a = c)] 
is the elevator rope acceleration. The last term in the 
righthand-side of 0 has been added to model the semi¬ 
active linear damper effect on the rope at the contact 
point 5(y — l + ldp), where <5 is the Dirac impulse 
function, kd p is the damping coefficient (the control 
variable) and Z^ p is the distance between the car-top 
and the point of attach of the damper to the controlled 
rope, e.g. [5]. 

The PDE 0 is associated with the following two 
boundary conditions: 

u{0,t) = fi{t) 

u(l(t),t) = f 2 (t ) W 

where /i(t) is the time varying disturbance acting on 
the rope at the level of the machine room, due to 
external disturbances, e.g. wind gust. f 2 (t) is the time 
varying disturbance acting at the level of the car, due 
to external disturbances. In this work we assume that 
the two boundary disturbances acting on the rope are 
related via the relation: 

/ 2 (f) = h(t)sin (*^ H ^ ), He R (3) 

where H is the height of the building. This expression 
is an approximation of the propagation of the boundary 
disturbance f\ along the building structure, based on 
the length l, it leads to f 2 = f\ for a length 0 (which 
















is expected), and a decreasing force along the building 
until is vanishes at l = H, f 2 = 0 (which makes sense, 
since the effect of any disturbance f\, for example 
wind gusts, is expected to vanish at the bottom of 
the building). As we mentioned earlier the tension of 
the rope T(y) depends on the type of the rope that 
we are dealing with. In the sequel, we concentrate on 
the main rope of the elevator, the remaining ropes arc 
modelled using the same steps by simply changing the 
rope tension expression. 

For the case of the main rope, the tension is given by 

T(y,t) = (m e + p(l(t)-y))(g-a(t)) + 0.5M cs g (4) 


where g is the standard gravity constant, m e ,M cs are 
the mass of the car and the compensating sheave, 
respectively. Next, we reduce the PDE model ([l} to 
a more tractable model for control, using a projection 
Galerkin method or assumed mode approach, e.g. [10], 
[ 11 ]. 

To be able to apply the assumed mode approach, 
let us first apply the following one-to-one change of 
coordinates to the equation (jT|) 

u(y,t) = w(y,t) + V fi(t) + (5) 

One can easily see that this change of coordinates 
implies trivial boundary conditions 


w( 0, t) = 0 
w(l(t),t ) = 0 


( 6 ) 


Now instead of dealing with the PDE ([!]) with non-zero 
boundary conditions, we can use the equivalent model, 
given by equation ([7]) associated with trivial boundary 
conditions Q. 

Following the assumed-modes technique, the solution 
of the equation ([7]), Q writes as 

j=N 

w (y,t) = X] qjiWji.y’t), a<en (io) 

j=i 


where N is the number of bases (modes), included in 
the discretization, (f>j. j = 1 ,..., N are the discretiza¬ 
tion bases and qj, j = 1,... A' arc the discretization 
coordinates. In order to simplify the analytic manipu¬ 
lation of the equations, the base functions are chosen 
to satisfy the following normalization constraints 


WW W(*) 

/ <p 2 j(y,t)dy = l, / <t> i (y,t)<t) j (y,t)dy = 0, Vi ± j 
Jo Jo 


( 11 ) 

To further simplify the base functions, we define the nor¬ 
malized variable, e.g. [5], [3] £(t) = , and the nor¬ 

malized base functions <f>j(y,t) = y , j = 1, ...,1V. 
In these new coordinates the normalization constraints write 

as /o^ 2 (£K = !» V * ^ 3 - After 

discretization of the PDE-based model ([7|, <[8]» and <0 (e.g. 
refer to [3]), we can write the reduced ODE-model based on 
A-modes as 


Mq + (C + CU)q + (K + ku)q = F(t) + F{t)U , 

where 


q £ R n 
( 12 ) 


After some algebraic and integral manipulations, the 
PDE model (jTJ) writes in the new coordinates as 

p & + MOpmjm. + {p y2 - T (v .*)) Sa 
+(G(t) + vk dp 6{x-l + l dp )) §f 

+(c p + kd p S(x — l + ldp))^f 
= y(-psi(t) - c p s 2 {t)) - pf[ 2) + s 4 (f) 

+ i;h(t))kdpS(x - l + l d p) 

-vkdp^^dtx - l + l dp ) 


(7) 


where G(t) = pa(t ) — +c p v{t), and the S{ variables 


are defined as 

Zi< 2 >-2Z 2 


S l(f) = ftVi=Af/l(f) + 2 T/ 1 
! (l 3 fi 2) -f2pP 2) +2li 2 h-2lHf 2 ) _ 

s 2 (f) = 4+ 4-/apt 


( 8 ) 


S3 (t) = h -J ± 

S 4 (f) = -2v(t)ps 2 (t) - G(t)s 3 (t ) - c p fi(t) 

associated with the two-point boundary conditions 

w(0,t) = 0, w(l(t), t) = 0. (9) 


U — kdp 

pr 2 i 2 6ij-pr 2 i 2 

+pl~ 1 {g-a{t)) fg( l-€)V>i ( 0 lp'j( 0 d £+Tn,,r 2 (g-a(t)) f* (f)V'j 

+ P q- 2 i 2 -z- 1 o(o.55 iJ -/ 0 1 (i-«)q i («)^(«)d«)- 

cjr^fo 1/H («)</>' (^idi+QXSiP+OXM^gl- 2 f* ^ (£)*/>'■ («)d? 

Fi(t)—-lVl(psi(t)+c p S2(t)) fg i/>i(£)(d(; 

+v / z(s 4 (t)-p/[ 2) (i)) foil>i(€)d£ 

A(t)=A^q^) + hd|£^(h7fii)(zz- 2 (/ 1 -/ 2 )+z- 1 (/ 2 -/ 1 )) 

r _ / 0 , i j- j 
°ij ~ 


1 , i = 3 


(13) 

where i,j £ (1,..., A}, and Sk, k = 1,2,3,4 are given in 

©• 

Remark 1: The model ( p~2] >, ( |T7| ) has been obtained for the 
general case of time-varying rope length Z(f), however, in this 
paper we only consider the case of stationary ropes l = cte, 
which is directly deduced from (12 1 , (13 1 , by setting 1 = 1 = 


0, Vi. The case where the car is static and the control system 























has to reject the ropes’ oscillations is of interest in practical 
setting. Indeed, besides the case of commercial buildings at 
night (where the elevators are not in use), in many situations 
where the building is swaying due to external strong weather 
conditions, the elevators are stopped, for the security of the 
passengers. The control system is then used to damp out the 
ropes sway, to avoid the ropes from damaging the elevator 
system, and be still functional after the external disturbances 
have passed. [[] 

III. Main result: Lyapunov-based semi-active 

DAMPER CONTROL 


The first controller deals with the case where the building, 
hosting a stationary elevator (stopped at a given floor), 
sustains a brief (impulse-like) external disturbance. For 
example, an earthquake impulse with a sufficient force to 
make the top of the building oscillate, or a strong wind 
gust that happens over a short period of time, exciting 
the building structure and implying residual vibrations 
of the building even after the wind gust interruption. In 
these cases, the elevator ropes will vibrate, starting from a 
non-zero initial conditions, due to the impulse-like external 
disturbances (i.e., happening over a short time interval), and 
this case correspond to the model ( 121 , (131 with non-zero 
initial conditions and zero external disturbances. We can 
now state the following theorem. 


Theorem 1: Consider the rope dynamics ( fl2) , ( fl3] >, with 
non-zero initial conditions, with no external disturbances, 
i.e., fi(t) = / 2 (f) = 0,Vt, then the feedback control 

a T Ca 

U(z) = Umax , _ = (14) 

y 1 + ( q T Cq ) 2 

where z = ( q T ,q T ) T , implies that q(t) —> 0 , for t —1 0 , 
furthermore \U\ < u max , Vi, and \U\ deceases with the 
decrease of q T Cq. 


Proof: We define the control Lyapunov function as 


V(z) = \q T (t)Mq(t) 


1 


q T (t)K{t)q(t) (15) 


where 2 = (q T , q T ) T . 

First we compute the derivative of the Lyapunov function 
along the dynamics (12 1 , without disturbances, i.e., F(t) = 
0, F(t) = 0Vi 

V(z) = q T (-Cq - CUq - Kq) + q T Kq 
= —q T Cq — q T CqU 1 ’ 


1 Of course, there are also practical cases were the cars are in 
motion and an external disturbance occurs. These cases correspond 
to a time-varying rope length, which we have also studied, however, 
due to space limitation, we could not include all the results in this 
paper. The case of time-varying rope length will be presented in 
another report. 


Next, using U defined in (14}, we have 


V(z) < —u 


max 


(i T C4) 3 

y/l + (qTCq ) 2 


( 17 ) 


Using LaSalle theorem, e.g. [12], and the fact that C is 
symmetric positive definite we can conclude that the states 
of the closed-loop dynamics converge to the set S = {z = 
(< q T ,q T ) T e R 27V , s.t. q = 0 }. Next, we analyze the 
closed-loop dynamics: Since the boundedness of V implies 
boundedness of q. q and by equation ( 121 , boundedness of 
q. Boundedness of q. q implies the uniform continuity of 
q, q , which again by © implies the uniform continuity 
of q. Next, since q —> 0, and using Barbalat’s Lemma, 
e.g. [12], we conclude that q —> 0 , and by invertibility 
of the stiffness matrix K + /3U we conclude that q —> 0. 
Finally, the fact that V is a radially unbounded function, 
ensures that the equilibrium point (q, q) = ( 0 , 0 ) is globally 
asymptotically stable. Furthermore the fact that \U\ < u max , 
and the decrease of \U\ as function of q T Cq is deduced from 

equation (14 1 . since . ^ c f= < 1 and . 1 . = < 1 . 

Yl + (9 T Cg) 2 - Yi + f'F’C?) 2 “ 


Remark 2: It is clear from equation © that the trivial 
choice of a constant positive damping control U, will also 
imply a convergence of the sway dynamics to zero. However, 
the controller ( [T4| has the advantage to require less control 
energy comparatively to a constant damping, since by con¬ 
struction of the control law ( fl4| >, the stabilizing damping 
force decreases together with the decrease of the sway. 

The controller U given by ( fl4| ) does not take into account 
the disturbance F(t) explicitly. Next, we present a controller 
which deals with the case of a static elevator in a building 
under sustained external disturbances, e.g. sustained wind 
forces on commercial buildings at night where the cars are 
static. In this case F(t) 7 ^ 0, F(t) 7 ^ 0 over a non-zero time 
interval, and satisfy the following assumption. 


Assumption 1: The time varying disturbance functions 
fi, /2 are such that, the functions F(t), F(t) are bounded, 

i.e. 3( Fmaxi Fmax') i S.t. \m\ — Fmax-i \m\ < 

F Vf 

x max ? v v. 

Theorem 2: Consider the rope dynamics ©, <©, un¬ 
der non-zero external disturbances, i.e., fi(t) 7 ^ 0 , / 2 (f) 7 ^ 0 
satisfying Assumption [T] with the feedback control 


U(z)—u . 
+q T Cq 
+q T Cq 


q C g 


yi+(i T Cj) 2 

v lmax {r'max I 9 I' 


,+fmoil!it 


y/l + (q T Cq) 2 (F ma 
v 2max (191 -^maa 

Vl + (9 T C<i) 2 (|<i|F„ 


max\q\u 

c lg|) 

n.ax l<?l) 2 


(18) 


where u maXpl v\max: U2max V* 0 , are chosen s.t. ( u m ax p 
Vi max T U 2 max) 5 u m ax and z = {c[ 7 q ) . Then, if we 


define the two invariant sets: 

Si = {(g T ,g T ) T eM 2W , s.t. 

iq T Cg) 2 

^Jl+(q T Cq) 2 (F ma . 1 : \q\v lrna .^+F Trla , ;c \q\u maa:p ) 2 

























Next, if we define the term 


, and 


S 2 = {(q T ,q T ) T e K 2A \ s.t. 

( fCq ) 2 

•\/l+(g T C , (j) 2 (|ij|F maa .+i;2 TOOa; -F , max |9|) 2 


< 


1 


} 


, the controller (18 i ensures that the state vector z converges 
to the the invariant set Si or S 2 , and that \U\ < u max , Vi. 
Proof: Let us consider again the Lyapunov function 
Its derivative along the dynamics ( fl2| ), with non-zero 
disturbance, i.e. F(t) 7 ^ 0, F 0, writes as 


V(x) = q T {-Cq - CUq - Kq + FU) 
+q T I<q + q T F(t) 

= ~q T Cq - q T CUq + q T FU + q T F(t) 
< - q T CUq + q T FU + q T F(t) 


Now, we use the concept of Lyapunov reconstruction, e.g. 
[13], we write the control 


U = U nomp +V 1 +V 2 (20) 

where, u nomp is the nominal controller, given by ( | 1 4[ i with 
Um ax = u maXp , designed for the case where F(t) = F(t) = 
0, Vf. The remaining terms Vi, v 2 are added to compensate 
for the effect of the disturbances F and F, respectively. We 
design v\ and V 2 in two steps: 

- First step: We assume that f / 0, F(t) = 0, Vf, and 
design vi to compensate for the effect of F. In this case 

U = U n0 m p T U\. 

In this case, the Lyapunov function derivative is bounded as 
as follows 

V{x) <-q T C(u nojnp +fi )q + q T F(u nomp -Mq) 

< ~q T Cviq + q r F(u nomp + i>i) 

< -q T Cqv 1 + q T Fu nomp + q T Fv 1 

which under Assumption [T] gives 

< ~q T Cqv 1 + \q\F max 'U'maXp “1“ \ ( l\F'max'V\ max 


Now, if we define (to simplify the notations), the term 


Ti — +\q\F ma x U'maXp “1“ \o\Fmax'V\ r, 


and if we choose the controller 

^ _ VlmaxTiq Cq 

iJl+T?(qTCq) 2 

This leads to the following Lyapunov function derivative 
upper bound 


<Ti(l 


VlmaUq T Cq ) 2 

\A +T?(q T Cq) 2 


- Second step: We assume that F(t ) 7 ^ 0, F(t ) f 0, and 
design v 2 to compensate for F. In this case we write the 
total control as U = u nomp +vi +V 2 - Now the upper-bound 
of the Lyapunov function derivative (along the total control) 
writes as 

< B\ — q T Cqv 2 + q T Fv 2 + q T F 

< Bi- q T Cqv 2 + \q\F max v 2 max T 


T 2 — \q\F m axU2max T 

and if we choose the controller 


V 2 


V2 max T 2 q T Cq 
V^l +Tf{q T CqY 


This leads to the following Lyapunov function derivative 
upper-bound 


<B 1 +T 2 ( 1 
< B\ + B 2 


V2 m aUq T Cq) 2 

\/l +T%(qTCqP 


where U 2 = F (1 
By choosing u-|. 


W2mo JjffCqf 


yJx+T^qrCqY 1 ' 

> u 2max high enough the two terms 
B\, B 2 will be made negative. We can then analyze two 
cases: 

1- First, the trajectories keep decreasing until they reach the 
invariant set 


51 = {(g T ,<f) T eK 2JV , s.t. Bi> 0} 

Then, the trajectories can either keep decreasing (if \B 2 \ > 
B\) until they enter the invariant set 

5 2 = {(q T ,q T ) T € R 2N , s.t. B 2 > 0 } 


or they get stuck at Si. 

2- Second, the trajectories decrease until they reach the 
invariant set S 2 first, and stay there, or keep decreasing until 
they reach the invariant set Si. 

In both cases, the trajectories will end up in either .Sj or S 2 . 
Finally, the fact that \U\ < u max is obtained by construction 
of the three terms, since from (14 1 , we can write \u nomp \ < 
Umax p and by construction \vf\ < Vi max , \v 2 \ < v 2 max, 
which leads to \U\ f Umax p F Vi max “h v 2 max S: Umax' ® 

Remark 3: We want to underline here the fact that, con¬ 
trary to the previous case of impulse disturbances, in this 
case of sustainable external disturbances, the use of a simple 
passive damper, i.e. U = cte, might actually destabilize the 
system. Indeed, by examining equation ( p~9| ), we can see that 
if the U is constant the positive term +q T FU, which is 
due the external disturbance, could overtake the damping 
negative term —q T CUq, leading to instability of the system. 

Remark 4: The controllers {I4]>, (T8]i are state feedbacks 
based on q. q, these states can be easily computed from 
the sway measurements at N given positions y( 1), ...,y(N), 
via equation (lOi. The sway w(y,t) can be measured by 
laser displacement sensors placed at the positions y(i), i = 
1,2 ,...N. along the rope, e.g.[14], subsequently q can be 
computed by simple algebraic inversion of ( fT 0 | ), and q can 
be obtained by direct numerical differentiation of q. 


















Parameters 

Definitions 

Values 

n 

Number of ropes 

8[ ] 

m e 

Mass of the car 

3500 [kg] 

p 

Main rope linear mass density 

2.11 [kg/m\ 

i 

Rope maximum length 

390 [to] 

H 

Building height 

402.8 [to] 

Cp 

Damping coefficient 

0.0315[At.sec/m] 


TABLE I 

Numerical values of the mechanical parameters 


IV. Numerical example 

In this section we present some numerical results obtained 
on the example presented in [1], The case of an elevator 
system with the mechanical parameters summarized on Table 
|T] has been considered for the tests presented hereafter. We 
write the controllers based on the model {12} , {TJ} with one 
mode, but we test them on a model with two modes. The 
fact is that one mode is enough since when comparing the 
solution of the PDE {7} to the discrete model ( p~2] > the higher 
modes shown to be negligible, and a discrete model with one 
mode showed a very good match with the PDE model, but to 
make the simulation tests more realistic we chose to test the 
controllers on a two modes model. Furthermore, to make 
the simulation tests more challenging we added a random 
white noise to the states fed back to the controller (equivalent 
to about ±1 cm of error on the rope sway measurement 
from which the states are computed, see Remark 1), and 
we filtered the control signal with a first order filter with 
a cut frequency of 10 hz and a delay term of 5 sampling 
times, to simulate actuator dynamics and delays due to 
signal transmission and computation time. First, to validate 
Theorem 1, we present the results obtained by applying 
the controller {14} , to the model {12} , {L3} , with non-zero 
initial conditions q( 0) = 20, g(0) = 5, and zero external 
disturbances, i.e. f±(t) = / 2 (i) = 0, Vi. In these first tests, 
to show the effect of the controller {14} alone, without the 
‘help’ of the system’s natural damping, we fix the damping 


coefficient to zero, i.e. c p = 0. We apply the controller (14 1 , 
with u m ax = 10 9 Nsec/m. Figures in show the rope 
sway obtained at half rope-length y = 195 m with and 
without control. Without control the rope sway reaches a 
maximum value of about 1.5 m. With control we see clearly 
the expected damping effect of the controller, which reduces 
the sway amplitude by half. The corresponding control force 
is depicted on Figures B0 We see that, as expected from 
the theoretical analysis of Theorem 1, the control force 
remains bounded by a maximum value of 40 kNsec/m, 
which is easily realizable by existing semi-active damper, 
e.g. magnotorheological damper. Furthermore, as proven in 


2 The figures' zoom is included for the reader to have a better 
idea about the signals shape. 
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Fig. 2. Rope sway at y — 195 m: No control (thin line)- With 
controller (|14[> (bold line) 



Fig. 3. Zoom of rope sway at y = 195 m: No control (thin line)- 
With controller {14} (bold line) 


Theorem 1, the control amplitude decreases with the decrease 
of the sway. 

Next, we consider the model {12}, {13} with no-zero 
disturbance signals: fi(t) = 0.2sin(2n.0.08t), and f 2 being 
deduced from /i via equation 0. We underline that we have 
purposely selected the disturbance frequency to be equal to 
the first resonance frequency of the rope, to simulate the 
‘worst-case scenario’. In this case we apply the controller 
(18 i, with the parameters u maXp = 10 9 Nsec/m , vi max = 
= 10 5 Nsec/m , F max = F max = 1. We show on 


^2 max 

Figures [6] [7] the sway signal in the uncontrolled and the 
controlled case. We see that the sway steady state maximum 
amplitude is reduced from 8.4 to in the uncontrolled case 
to 2.4 to with control. The noisy (due to the simulated 
measurements noise) bounded and continuous control signals 
are reported on Figures [8] [9] 
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Fig. 4. Output of controller (|14[) 



Fig. 5. Output of controller 03 - Zoom 



Time [sec] 


Fig. 6. Rope sway at y = 195 m: No control (thin line)- With 
controller (fl8l) (bold line) 



Fig. 7. Zoom of rope sway at y = 195 m: No control (thin line)- 
With controller {18} (bold line) 


V. Conclusion 

In this paper we have studied the problem of semi active 
control of elevator rope sway dynamics occurring due to 
external force disturbances acting on the elevator system. 
We have considered the case of a static car, i.e. constant rope 
length and have proposed two nonlinear controllers based on 
Lyapunov theory. We have presented the stability analysis of 
these controllers and shown their efficiency on a numerical 
example. The semi-active stabilization problems related to 
time-varying rope lengths, i.e. moving car, will be presented 
in a future report. 
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Fig. 9. Output of controller l|18|l- Zoom 
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